Tunneling spectroscopy using a probe qubit 
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We describe a quantum tunneling spectroscopy technique that requires only low bandwidth con- 
trol. The method involves coupling a probe qubit to the system under study to create a localized 
probe state. The energy of the probe state is then scanned with respect to the unperturbed energy 
levels of the probed system. Incoherent tunneling transitions that flip the state of the probe qubit 
occur when the energy bias of the probe is close to an eigenenergy of the probed system. Monitoring 
these transitions allows the reconstruction of the probed system eigenspectrum. We demonstrate 
this method on an rf SQUID flux qubit. 
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Recent technological advances have allowed the con- 
struction of mesoscale systems of individual quantum 
elements, including hundreds of trapped ionsi, 14 en- 
tangled ions^, nanomagnetic systems assembled out of 
magnetic atoms on metallic surfaces^, ultracold *^Rb 
atoms in optical lattices^, and arrays of superconducting 
devicesS,jl. While the study of small numbers of atoms or 
devices often involves direct manipulation and full state 
tomography, these techniques become impractical in the 
mesoscale regime. As a result, there is a need for tools 
that are applicable when one has a mesoscale system with 
limited control over its individual elements. 

Tunneling spectroscopy is a powerful tool for study- 
ing condensed matter systems. It has been used to 
push the limits of our understanding of many-body 
physics, as in recent studies of two dimensional elec- 
tron systems in high magnetic field using time do- 
main capacitance spectroscopy^ and scanning tunneling 
spectroscopy^. Tunneling spectroscopy can also be used 
to directly validate numerical or analytical models of 
complex systems, such as the single particle states of 
CdSe quantum dots^ or the electronic wavefunctions of 
carbon nanotubes^. Motivated by the ability of tun- 
neling spectroscopy to probe the quantum behaviour 
of mesoscale systems, we have developed an analogous 
method that is applicable when one has limited control 
over a large system. The large system is probed using 
a dedicated probe qubit with its own readout and low 
bandwidth control of its Hamiltonian. We have termed 
this new technique qubit tunneling spectroscopy (QTS). 

A related method has been proposed in Ref. [ll| where 
the probe qubit must be perturbatively coupled to the 
system under study. In QTS, the requirement for this 
weak coupling has been removed through the use of a 
compensation bias (as explained below). Further, the 
algorithm of Ref. 11 is designed to operate on a gate 
model quantum computer while we demonstrate QTS on 
a system with much more limited control. 

QTS requires a probe qubit P that can be described 



by a generic two-level system Hamiltonian: 

Hp = --Apo-x^P - -ep&z^p, (1) 

where a^^p and ct^^p are Pauli matrices operating on P, 
and both parameters ep and Ap should be controllable. 
The eigenstates of az,p with eigenvalues -1-1 and -1 are 
|t)p and |i)p, respectively. The <Tz,p eigenstates should 
be distinguishable by a readout mechanism. Let there be 
a system S, governed by some Hamiltonian Hs, that one 
would like to study. 

To perform QTS we require a coupling between the 
probe qubit and a parameter of the system (described 
by an operator C) as well as a controllable compensation 
bias ecomp coupled to that same parameter. In this case, 
the system plus probe Hamiltonian can be expressed as: 

Hs+p = Hs + Hp + Jcjz.pC + -ecompC*, (2) 

with J the strength of the probe qubit-system interac- 
tion. 

For general ecomp the eigenstates oi S + P are not rep- 
resentative of those of S. However, in the special case 
Ecomp = —2 J, the spectrum of Hs+p splits into two 
qualitatively different manifolds, M-f and Ai^, wherein 
the probe qubit P is in state |t)p and \i)p, respectively. 
Hs+p can then be rewritten as: 

Hs+p= (^epIs-2JC + Hs)(E>\i)p{i\p (3) 

+ ^S® lt)p(tlp 

- Y^5® (i;)p(tip + it)paip) 

where Is is the identity operator on system S. In the 
case where the third line of Eqn.Ois perturbatively small, 
the first line is the Hamiltonian of Ai^ and the second 
line that of M^. Thus, the energy spectrum of is 
identical to that of Hs. Further, the first term of the 
first line of Eqn. [3] shows that the energy of all states 



2 



IT)Qli>^ 



i4J 



7;(IT)g-|i)g)IT)p 



\1)q\1)p 
Ml 



Toi ^(IT)g + li>g)IT>p 



FIG. 1. Energy level diagram for the Q + P qubit system 
described by the Hamiltonian in Eqn.|4]at eq — and ecomp ~ 
—2J. The system separates into two distinct manifolds, 
and Mf, that differ in the orientation of the probe qubit 
P. Spinor notation indicates state of qubits Q and P on 
the left and right, respectively. Allowing for weak tunneling 
in qubit P facilitates incoherent tunneling processes Foi and 
ro2 between the manifolds. QTS proceeds by initializing the 
system in the lowest energy state of 7V(4.and finding values of 
e p for which transition rate F peaks due to resonant processes 
Foi and F02. 



in Ml can be shifted with respect to those of A^^ by 
adjusting the probe energy bias ep. For small enough 
Ap, the third line of Eqn.[3]gives rise to incoherent inter- 
manifold tunneling^^ between any state \k'^) |4-)p of Aii 
and any state \k) |t)p of A^^ with a rate proportional to 
|Ap(fcl|fc>P- 

The QTS method begins by initializing the system 
into the lowest energy state of A^j^. The tunneling rate 
between manifolds peaks when an eigenstate of M.-f is 
brought into resonance with the initial system state in 
A^^ by adjusting ep. This resonant tunneling transi- 
tion between manifolds flips the state of the probe qubit, 
which can be easily detected. Thus, to perform QTS, 
one measures the initial transition rate F = \dP/dt\^^f^, 
where P is the probability of observing the probe qubit 
in its initial state, as a function of ep, the probe energy 
bias. Scanning ep and locating peaks in F allows one to 
map out the eigenspectrum of which is identical to 

that of Hs if the compensation bias ecomp is set to —2 J. 
Note that errors in this compensation bias will skew the 
energy spectrum of the probed system. The errors in the 
extracted energy spacings of the spectrum are bounded 
by the compensation bias error. 

To experimentally demonstrate QTS we take a target 
system comprising a single qubit Q, governed by a Hamil- 
tonian Hq (Eqn. [1] with P — > Q). We couple the probe 
qubit P with strength J through a mutual CTz interaction 
to qubit Q. In this case, Hs Hq and C — > a^^g in 
Eqn. [21 yielding 
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Q+P 



Ho + Hf 



(4) 



Ap is chosen to be small compared with all other terms 
so the eigenstates of the probe are to good approximation 
|t)p and \\) p. Setting tcomp — — 2J then yields a Hamil- 
tonian of the form given in Eqn. [3] For the particular case 



eg = 0, we show how the theoretical energy spectrum of 
the coupled two qubit system splits into two manifolds 
in Fig.[TJ Note that the eigenstates of shown therein 
are superpositions of the W)q and \\.)q states. Changing 
ep allows the lowest state in A^^, |J,)p, to be brought 
into resonance with the states in M.^j. The system plus 
probe can tunnel from M.^ to M.^ through the incoherent 
processes labelled as Fqi and F02. 

The two qubits, probe P and target Q, used in the ex- 
periment were rf SQUID flux qubits on a D-Wave quan- 
tum annealing processor thermalized to a temperature 
T = 12 mK. A description of a chip similar to that used 
in this study can be found in Ref. Il3l . The low energy rf 
SQUID flux qubit Hamiltoniarti^ has a direct mapping 
onto Eqn. [T] 
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where we have performed the substitution Aq 



and eg — 2 



$5,, with $ 



and 



<i?Q being externally controlled flux biases and 



being the magnitude of the qubit persistent current. Note 
that both Ag and eg (through Iq ) are functions of 

<j,«^jj . The functional forms of these dependencies are 
determined by the physical parameters of the rf SQUID, 
as described in detail in Ref.llJ. If one considers the qubit 
Q as an Ising spin, then Ag corresponds to a transverse 
magnetic field, Iq is the magnitude of the spin, and $q 
is an applied longitudinal magnetic field. The physical 
Hamiltonian for the probe qubit Hp is found by replac- 
ing Q by P in Eqn. [5l The probe qubit had a persistent 
current \Pp\ = 1.0 ^A and Ap//i - 1 MHz. The small 
Ap was chosen so that the transition rate F of the probe 
qubit was contained within the dc to 3 MHz bandwidth 
of our slow control lines and to satisfy the incoherent 
inter-manifold tunneling condition. An on-chip tunable 
coupler— between the two qubits was programmed to at- 
tain an interqubit mutual inductance M = 2.0 pH. The 
resulting form for J in Eqn. |3]is J = M|/q||/p|. With 
these parameters, the spectral gap in the A^j, manifold, 
as depicted in Fig.[Tl satisfied 4J = 4M|/^||/];,| > ksT 

over the range of encountered in these experiments. 
Consequently, there was negligible thermal activation out 
of the initial state |i)g|4,)p to higher levels within A^j^. 
With these parameters, the compensation bias is explic- 
itly ^ecompO-z,Q = -J^z,Q = -M| 1 1 J|, |(T^,g and is ap- 
plied by adding an offset A/ 1 /p | ~ 1 m<I>o to the flux bias 
$Q of qubit Q. This compensation bias requires only 
careful calibration of probe parameters. 

The experimental method for initialization and read- 
out is the same as the two-qubit cotunneling technique 
described in Ref. [1^ The experiments described herein 
differed from the cotunneling experiment in three re- 
gards: First, in QTS one intentionally sets Ap ^ Ag, 
thus exploring an extreme limit of the mismatched tun- 
neling energy configuration described in Ref. [l5l . Second, 
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FIG. 2. (a) QTS tunneling rate data versus ep for <1>q = 
0.1m<l>o and $q"'"' = 0.637. Peaks in F are readily asso- 
ciated with the processes Fqi and denoted in Fig. [1] 
The centers of the peaks are found by fitting to a sum 
of two Gaussian peaks, (b) Qubit Q energy spectra ob- 
tained by QTS at control bias values: <1)q^^7$o = 0.640; (c) 
$Q-'-'/"I>o = 0.637. In both plots, the ordinate is the probe 

bias energy ep = 2|jp|"I>p and the abscissa is the flux bias 
$Q oc tQ applied to qubit Q. The grayscale indicates the 
transition rate F in the (3>Q,ep) plane. White circles denote 
the centers of peaks in F found by the Gaussian fits. Avoided 
crossings between two localized states, explicitly labelled in 
(b) as |t)(5 and \\.)q, are visible, (d) The difference in probe 
energy Sep between the two peak centers as a function of <l?^ 
for the datasets in (a) and (b). Results have been fit to the 
dispersion of the Hamiltonian in Eqn. [5] using Ag and | Iq | 
as free parameters. 



in QTS we use relatively large offset biases Ccomp/ j-fgl = 
lni$o in order to satisfy the compensation condition em- 
bodied in Eqn. [31 Third, whereas the dynamics studied 
in Ref. [T5| involved incoherent tunneling of the pair of 
qubits between localized initial and final spin states, in 
QTS the final state can place qubit Q in a delocalized 
(superposition) state, as depicted in Fig. [TJ 

A scan of the initial transition rate F versus ep at 



Q 



and <^''^" = 



0.637 is shown in Fig. [Sfa). The 
data clearly show two distinct peaks. These peaks are 
readily identified as the processes Fqi and F02 indicated 
in Fig. [TJ We fit such scans to a model composed of a 
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FIG. 3. Pq{^q'') vs l^Qi^^"^") from the Q spectra for a 
range of settings of the control bias $q''''. The data points 
from the fits in Fig. [2](d) are labelled with arrows. The data 
also have horizontal error bars approximately the size of the 
symbols. Results have been fit to the rf SQUID model from 
Ref. [H (solid line) with the rf SQUID capacitance C as the 
only free parameter. 



pair of Gaussian peaks in order to locate the peak cen- 
ters. Example maps of the initial transition rate F versus 
ep for a range of qubit Q flux biases around $q = are 
shown in Figs. [H^b) and (c) for two values of the control 
fiux bias $q"'"' (and therefore Ag). For clarity the cen- 
ters of the Gaussian peaks have been indicated by white 
circles. The peak positions reveal the avoided crossing 
between localized spin states W)q and ||)g. In Fig.[2Ud), 
we summarize the difference in probe energy 5ep between 
the two peaks as a function of $q . We then fit those re- 
sults to the eigenspectrum of the Hamiltonian in Eqn. [SJ 
using Aq and |/g| as free parameters. 

In order to crosscheck our QTS results, we have 
repeated the measurements and analysis that led to 
Fig. [Hd) for several values of ^q'^'' ■ By doing so, we 



generated maps of the qubit parameters Ag and 



I'r 



a function of 



A plot of the relationship between 



and Ag is shown in Fig. [3] This curve is completely 

determined by the rf SQUID inductance L and capaci- 
tance C of qubit Q. We have fit these results (solid line 
in Fig. [31 to a physical rf SQUID Hamiltonian (Eqn. 4 
in Ref. [ij) taking L = 355.5 pH, as determined by inde- 
pendent measurements, and using C as a free parameter. 
The best fit returned C = 118 ± 2 fF, which is a physi- 
cally reasonable value, given Josephson junction sizes and 
qubit wiring geometry. The single parameter fit models 
the data well, implying that QTS has correctly extracted 
the low energy spectrum of rf SQUID qubit Q. 

QTS could, in principle, yield more than just the eigen- 
spectrum of system S. In particular, there is signifi- 
cant information contained in the spectral weight of the 
peaks in F. For example, the spectral lines inferred 



4 



from Fig. [IJb) are less pronounced for the upper level 
at $Q < and for the lower level at 4>q > 0. This is due 
to the proportionality of the initial transition rate F to 
the small overlap of the initial state of Q IDq ) with its 
final state (~ |t)Q)- Choosing ecomp = +2 J, instead of 
—2 J, yields a system in which the states in Ai^ exchange 
roles, thus yielding a new initial state \t)Q p- Repeat- 
ing the QTS experiment with this configuration should 
then swap the regions of high and low peak visibility seen 
in Fig. [2ljb). Thus the spectral weight contains informa- 
tion about the wavefunction of the probed system. 

Further information could be gleaned from the line- 
shapes of inter-manifold tunneling processes. We chose 
to fit tunneling rate peaks to Gaussians as we had an- 
ticipated that their lineshapes would be dominated by 
the incoherent tunneling of the slow probe qubit P, as in 
Ref. 16. A detailed analysis of the physical mechanisms 
that lead to particular lineshapes is currently underway. 

We have demonstrated a low bandwidth method, 
termed qubit tunneling spectroscopy (QTS), by probing 
the energy spectrum of a first qubit by using a second 
probe qubit to split the two qubit system into two man- 



ifolds of qualitatively different states. Transitions be- 
tween the two manifolds are monitored as a function of 
the energy bias of one of the manifolds. Transition rate 
peaks correspond to the presence of eigenstates in the 
target manifold under study at that energy bias. We 
validated this method by verifying that rf SQUID flux 
qubit energy spectra measured in this manner are con- 
sistent with an rf SQUID Hamiltonian. QTS provided a 
direct measurement of the first qubit 's energy splitting 
Aq that was three orders of magnitude larger than the 
measurement bandwidth Ap. 

While the demonstration in this paper was limited to 
a single qubit, QTS is extensible to larger numbers of 
qubits and to other physical systems, provided one has 
good control and readout of the probe and a method of 
applying a compensation bias. We anticipate that QTS 
will be a valuable tool for studying mesoscale systems. 
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